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Настоящая рукопись является первой из публикаций, посвященных динамике вытеснения нефти 
водной эмульсией. Предполагается, что обе жидкости разделены неизвестной поверхностью и 
их динамика описывается системой уравнений Стокса для вязких сжимаемых жидкостей. Одна 
из жидкостей (нефть) является каплей, окруженной другой жидкостью (водной эмульсией). На 
границе раздела двух жидкостей выполнены стандартные условия непрерывности перемещений 
и напряжений и дополнительное условие того, что область занятая каждой из жидкостей яв- 
ляется индивидуальным объемом, то есть состоит из одних и тех же частиц. Следуя идеям О. 
Ладыженской, Н. Уральцевой и В. Солонникова исходная задача разбивается на ряд последова- 
тельно решаемых задач, первой из которых является модельная задача дифракции для системы 
Стокса для двух полупространств. В ней необходимо получить решение в явном виде некоторые 
сингулярных интегралов. Далее для произвольного интервала времени [0, Т] решатся задача 
дифракции, когда форма капли есть заданная гладкая замкнутая ограниченная поверхность и, 
наконец, решается исходная нелинейная задача об определении формы капли в процессе ее дви- 
жения. 


Ключевые слова и фразы: Задачи со свободными границами, задача дифракции для системы 
уравнений Стокса, представление решений в канонических областях в форме сингулярных инте- 
гралов, метод Ньютона-Канторовича. 


ПОГ: 10.31857/52686954722020084 


ВВЕДЕНИЕ И ФОРМУЛИРОВКА РЕЗУЛЬТАТОВ 


Задача о вытеснении нефти водной эмульсией в горных породах является основной при создании 
гидродинамического симулятора нефтяного месторождения. Естественным является предположение 
о том, что соответствующая математическая модель должна базироваться на постулатах классиче- 
ской механики Ньютона сплошных сред [1], [2]. С другой стороны, во всех существующих гидро- 
динамических симуляторах нефтяных месторождений ("Эклипс"и "Наск ОП"фирмы ЭсШашфегеег, 
"Тешрез"фирмы "Вохаг "УТР"фирмы Гапашагк и "ТипейУ"фирмы "ЭбапЧага ОП ап Тгизё") базо- 
вой математической моделью является математическая модель Баклея-Леверетта [3], которая, по своей 
сути, является набором постулатов, никак не связанных с классической механикой Ньютона. Следуя 
В. Вигиаее, .. В. КеПег [4] и Е. Запсре?-Раепсла [5], для точного математического моделирования 
физических процессов в горных породах в первую очередь необходимо описать физический процесс 
законами классической механики Ньютона на микроскопическом уровне (характерный размер- десят- 
ки микрон), которые проверены опытом в течении столетий. Но такое точное описание совершенно 
бесполезно с практической точки зрения, поскольку любые численные реализации точных моделей 
потребуют месяцы и месяцы вычислений. Поэтому следующим шагом в методе этих авторов стал про- 
цесс усреднения (гомогенизации), позволяющий на основе соответствующей микроскопической модели 
точно описать физический процесс уже на макроскопическом уровне (характерный размер- десятки 
сантиметров или метров). Метод усреднения позволил построить точные модели для многих фи- 
зических процессов, когда микроструктура сплошной среды заранее известна. Но в случае задачи 
вытеснения нефти суспензией никакого прогресса не произошло. Во-первых, не было соответствую- 
щей математической модели на микроскопическом уровне и, во-вторых, если она бы и была, не было 
методов решения задач с неизвестными (свободными) границами и, тем более, методов усреднения 
таких задач. Задачи со свободными границами являются одними из наиболее трудных задач в теории 
дифференциальных уравнений с частными производными [6], [7]. Как мы уже отмечали выше, нанта, 
основная задача, разбивается на ряд последовательно решаемых задам: 
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[-модельная задача дифракции для уравнений Стокса в полупространствах; 

П-доказательство существования классического решения задачи дифракции для уравнений Стокса, 
в областях с заданной гладкой границей 5 раздела областей ®,, содержащей нефть, и ®.,, содержащей 
суспензию. 

Ш-доказательство существования классического решения движения капли нефти в водной суспен- 
зии в малом по времени (метод Ньютона-Канторовича). 

ТУ-доказательство существования классического решения движения капли нефти в водной сус- 
пензии в целом по времени (сведение исходной задачи к канонической задаче в полупространствах с 
помощью метода декомпозиции). 

У-усреднение. 

В настоящей публикации мы ограничимся этапами Г — ГУ. 

Исходная нелинейная задача дифракции для искомой (свободной) границы Г({) раздела двух вязких 
жидкостей состоит в нахождении вектора скорости %о, вектора перемещений о и давления ро нефти 
и вектора скорости %., вектора перемещений 12. и давления р. суспензии, и самой границы Г(® по 
следующим уравнениям 


У. (ам. В (т,9;) - 71) =0, (1) 
1 Эр; МИ 
<, т +У. 9; = 0 (2) 


со 
в областях Оз; = (—] О; (Ю, 7 = о, 8, граничным условиям 


#0 
(и В(т,5ъ) — ро!) < п >= (и, В(т, 5.) -р.П <п>, (3) 

дм; | 
о = 1, -эв = 20,8 (4) 


на искомой границе Г(#) и начальным условиям 


0; (2;0) = 122), ру(2;0) = 00(=), 1 =0, в, Г(0) = Го. (5) 


Очевидно, что траектории 1 = Х (&, +) транспортного уравнения 


Хх 
Ч = ху, Х (&;0) = &, ё=Е9. 0), 1=0, 5, (6) 


которые начинаются в точке & е 9 ;;(0), будут находиться @;;(#), 7 = 0, $ для всех # > 0. 
То есть, мы можем определить Г({) как поверхность 


Г) = {26 0:2= Хо(6), ЕТ}, 3 =ъ(Хозв, Хо: 0) =& ЕЕ Г. (") 


1 
В (1) - (7) О(т, и) = 5 (Уи + У*и) есть симметрический градиент вектор-функции и, [ единичный 
тензор второго порядка, ау и ср;, 1 = 0, $ безразмерные вязкости и безразмерные скорости звука 


соответственно в нефти и суспензии. 


Введем следующие обозначения [8]: 
1 1 | у? 5 — 92|? 
п 


(т; у’) = 
2 ) 4 [теб и |2 + у 


-функция Грина оператора Лапласа, в пространстве ®? и 


1 2% 1 
Сз(ж;у) = 1 (9 + 3 — 3?) *- 1’ |+ [3 +31)? (9) 


-функция Грина оператора Лапласа, в пространстве 3. 
По построению функции Грина бесконечно дифференцируемы при х = уи 


Сь(ж; у) = Сь(у; =), У„Сь (аж; у) = У,Сь(ж; у), К=12. 
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Кроме того положим 


15(и)(=) = | С2(жзу)щу У, #' Е В®, 
В2 


в(0е) = ||, Сы оби), 


(Ра) = ] ‚быв у) Кубы, зу В, 0, (10) 


КЗ — {т — (та, 2, хз), | Би (212+ |221? + [53]? < со}, КЗ — {т Е ВЗ, 3 > 0}, КЗ = {т Е ЕЗ, тз < 0}, 
х' = (21, 2) Е 2, |х'| = \/21 +22. 
Теорема 1. Функции 


шь() = ев) +..(5)@), +688, (11) 
и (#) = эВ.) +..0.)@), «ев (12) 


являются решениями задачи Дирихле для уравнения Пуассона 
Ди; = }, (2), жЕ Е, м”, 0) = (=), =, 8 (13) 


в полупространствах В 3, } = 0, 3. 
В частности, 
д / / / д / / / 
—— 13(05)(%',0) = -ио(ж',0) = -0.(%), —13(0,)(%',0) = и.(%,0) = 0, (2). (14) 
д тз дтз 
Если 1; Е С. (3) ц О; Е С° (В?) тощз Е С2+® (85), 1 = 0, з [1]. 
Мы используем обозначения функциональных пространств и норм в них, принятых в [10] ив [11]. 
Через С° (Я), & = 0, 1,2, ... мы обозначаем пространство С^*® (9) П 1..ъ (Я). 
Теорема 2. [8| Функция 


ща’) = ] ба") Пиау = 5) (15) 

Е2 

является, решением уравнения, Пуассона 
Д’и = {1 (2), х' Е В? (16) 


во всем пространстве В?. 
Если ЁЕ С“ (В?), тоие С2+° (В?) [10]. 
Теорема 3. Функции 


и (=) = —1з(Но) (2) + Ло (4о) (2) (17) 


из (2) = 13(Н+) (2) + 7. (7+) (2) (18) 


являются решениями задачи Неймана для уравнения, Пуассона 
Ди; = //(=), Е 3, 2 (2,0) = НД’), }=0, 8. (19) 
ТЗ 


Если }; Е с“ (В?), то и; Е С° (3), 7 = 0, 8 на произвольном интервале времелии [0, Т] [10]. 
Теорема 4. Задача (1) - (6) имеет, единственное классическое решение ‘шу, и; Е С31“(В3) иру Е 
С}, ] = 0, 3 на произвольном интервале времени [0, Т]. 


Теорема 5. Задача (1) - (7) имеет единственное классическое решение шу, и; Е С21“(9 (1) и 


РЕ С° (0, ;(1)), 1 = 0, 3 на произвольном интервале времени [0, Т]. 
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ДОКАЗАТЕЛЬСТВО РЕЗУЛЬТАТОВ 


Доказательство теорем 1-3 достаточно сложное и требует подробных вычислений, что не позволяет 
формат нашей статьи. Поэтому ограничимся идеями доказательств. Основные идеи доказательства 
состоят в замене краевого условия (7) на свободной границе транспортными уравнениями 

аа, _ данл 
и) — 9 Чи, — 0) — 5 
и Е Убив =0, а (0 = р 7=0, 5 (20) 
решение задачи дифракции в полупространствах и представление решения в виде сингулярных инте- 
гралов, использовании метода Ньютона-Канторовича |9] для доказательства существования классиче- 


ского решения задачи в малом по времени [9] и применение принципа декомпозиции, предложенного в 
работах [10] и [11] О. А. Ладыженской, Н. Н. Уральцевой и В. А. Солонникова, для исследования диф- 
ференциальных свойств в целом по времени классических решений задачи дифракции для системы 
уравнений Стокса в заданных областях. 

А именно, с помощью разбиения единицы исходная задача сводится к решению конечного числа, 
задач в малых областях с центром в точках жк Е 7%. Далее для каждой малой области участок 
границы выпрямляется переходом к новым локальным координатам у и задача в каждой малой области 
сводится к канонической задаче дифракции (1) - (6) с заданной правой частью }; Е С%. (ЕЗ) для новых 
искомых функций и,(у;%) = %;(щ;8), 9; (у; = р,(®;№ в полупространствах ие с границей раздела, 
Г = {из =0}: 


Решение последней дается формулами 


д д 
Ио ВИ) и = ВИ) +, 
к) + вк) + 


1 
Оз(уз 8 = 543( 


9 $оз й 0 $3 
О уз О уз 


Эт, Е 
а в ( —_ АЕ) + о )), (> — 13 ( — А! 13(02) + т И 


(у; 5, 


2 ‘ду дз’ 
Ко Оо 9 фо 3 Из 951 9 +4 3 
3(01) = = 15(03) + 1 — - - 1. а 
(071) Эл 2(Из) + и. О би, и + и. 25 дл 
[9] Ко Оо.2 О1оз Из 91.2 9 1з.з 
3(02) = 15 (Оз) + Г. Е | 1. ее —), (21 
3 2) Оу 2 3) Ко + Из 2 бу ду2 Ко + Из - уз ду2 
2 
0 СЯ 
а; (у;8) = 4; (у) ехр(-- 7, 1=0, 5. (22) 
й 
Здесь Ок(у'; 1) = изк(у’, 0; 8) = ок (у, 0;%), {; = УИХ,;), Л =0, зи \(а',0;1) = $,(2',0;#) = $(2'; 8). 
Обращаясь к [10] заключаем, что 
2-+а < а 
нрах, [М вз” < М ах, [Рева (23) 
где постоянная М не зависит от Ти {. 
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пбегРасе. ТБе зоПа зк@ебоп Баз а гуеп рено с этисбаге у а Аппепзющезз роте з17е = < 1. 6 за 
ЁВее БочпЧагу ргоешт, зшсе ш %Ъе писгозсор1с 4езсирйоп $Ъе пфетЁасе Бебмееп о| ап4 зизрепз1юй плиз Бе 
Чебегите4. Зась ргоЫетаз ате алпопя Ве 110$ ЧИ си ргоетаз ш Ве Ъеогу оЁ ратНа] Ч1егепйа] едиа 013, 
апа, аз а ге, ех1з6епсе гезиз аге роз е ошу 1осаПу ш ие. У\е Ш оба 21оБа] ш те гези ® Бу гедасше 
{Ве Ёее Бопп4агу ргоеш фо Ве ргоешт оЁ ВпАте %Ъе у1зсозКу оЁ Палаз, хасЬ у Ш Бе 4езстЬеа Бу %Ъе 
{тапзрог$ едаа1от. ТБе тала ригрозе оё $13 атае 1$ фо 4езсге %Пе }о1ё шовоп оЁ а зше о! агор ш $Ъе 
зиггоипА те уафег зазрепз!оп. 


Кеушотаз: ее Бочи4ату ргоетз, @ЯтасЯоп ргоетаз Юг ЭфоКез едиа оз, гергезетба оп оЁ зоаюопз Ш 
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